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Chiral transition and mesonic excitations for quarks with thermal masses 
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We study the effect of a thermal quark mass, mr, on the chiral phase transition and mesonic 
excitations in the light quark sector at finite temperature in a simple chirally-symmetric model. We 
show that while nonzero mr lowers the chiral condensate, the chiral transition remains of second 
order. It is argued that the mesonic excitations have large decay rate at energies below 2mr, owing 
to the Landau damping of the quarks and the van Hove singularities of the collective modes. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) at nonzero temper- 
ature exhibits rich physics, and is expected to undergo 
chiral and deconfinement transitions at a critical tem- 
perature T c . The RHIC experiments revealed interesting 
features of a quark-gluon plasma (QGP) phase above T c 
[l| . One of the historical problems is to explore the ex- 
istence and the nature of the mesonic excitations in the 
QGP phase 0, d, 0, [||. Recent lattice QCD simulations 
suggest the existence of such excitations in the heavy 
quark sector (2j. In the light quark sector, it is theo- 
retically predicted that there exist light collective modes 
in the sigma- and pi-channels, as the soft modes of the 
chiral phase transition Q. 

It is also known that the quarks show a peculiar be- 
havior in the QGP phase. At extremely high temper- 
atures, the hard thermal-loop (HTL) approximation Q 
can be adopted, and the quark propagator has two col- 
lective excitations, normal quasiparticles and plasmino 
modes [1, 0, H|. The excitation spectra of both collec- 
tive modes have a mass gap termed the "thermal mass" 
rriT ~ gT, where g and T are the gauge coupling and 
the temperature, respectively. Unlike the Dirac mass, 
the thermal mass does not break the chiral symmetry. 
Although the validity of the HTL approximation may be 
violated around T c , lattice simulations suggest a large 
quark mass near but above T c [9(. If the quarks have 
large thermal masses near T c , they can affect the proper- 
ties of mesonic excitations, especially in the light quark 
sector Q , and also the chiral transition at finite temper- 
ature. We notice that the quark spectrum can have a 
thermal mass in the vicinity of T c as a result of the soft 
mode of the chiral transition [l(| ■ 

In the present work, we explore the effects of the ther- 
mal mass on the chiral phase transition and mesonic exci- 
tations at finite temperature in a simple model. We em- 
ploy the propagator obtained in the HTL approximation 
for the quarks, with the thermal mass my introduced as 



a parameter. We also adopt a chirally-symmetric four- 
Fermi interaction in the scalar and pseudoscalar chan- 
nels. Using this model, we show that the order of the 
chiral transition does not change, while the value of the 
chiral condensate is suppressed, when tot ^ is in- 
cluded. It will be shown that even at finite tot, there 
are soft modes of the chiral transition in the scalar and 
pseudoscalar channels [||]. Using these modes, we study 
the effect of tot on the mesonic excitations above T c . 

How does the thermal mass of the quark affect the 
properties of mesonic excitations? Naively, a mesonic 
excitation below the expected threshold energy, w t hr = 
2tot, is tightly bounded, since the decay process into 
quark and anti-quark is forbidden below Wthr- We show 
in the present work, however, that the effect of the ther- 
mal mass is completely contrary to this statement; if the 
quark has the spectrum calculated in the HTL approx- 
imation with a thermal mass tot, the decay rate of the 
mesonic excitations is enhanced below the energy 2tot- 
It is discussed that these decay rates mainly originate 
from the continuum in the quark spectrum which phys- 
ically corresponds to the Landau damping of the quasi- 
quarks. 

This paper is organized as follows. In the next Section, 
we present the model employed in this work. In Section 
III, we calculate the chiral condensate in the mean-field 
approximation and demonstrate the chiral phase transi- 
tion in our model. We then calculate the spectral func- 
tion in the scalar and pseudoscalar channels above T c 
in Section IV. Section V is devoted to discussions and a 
summary. 



II. MODEL 

Before we present the model employed, let us first 
briefly review the property of the quark propagator in the 
high temperature limit of QCD, in order to see the emer- 
gence of the gapped collective excitations in the quark 
spectra. In this limit, the retarded quark propagator is 
calculated by the HTL approximation as 



'Electronic address: 
t Electronic address: 



hidaka@quark.phy.bnl.gov 



-fHTL , 



[(cj + 177)7° - p • 7 



S 



HTL/ 



if], p)} 



kitazawa@quark.phy.bnl.gov 



(1) 



2 



where p = \p\ and 



iHTL 



(u,p) = ^Q - 7° + — 1 - -Qo - 7 • P, 
P \pj P P \pj 

(2) 



is the quark self-energy in the one-loop level with m\ = 
(l/6)g 2 T 2 and Q = (1/2) ln(x + l)/{x - 1) @]. Notice 
that the quark propagator S 111 " 1 " (uj, p) is chirally symmet- 
ric as one can easily check that it anti-commutes with 7 s . 
The quark propagator Eq. (fTJ) can be decomposed as 



S^(u,p) = (S^(oj,p)P + (p) + S^(c,p)P_(p)) 7 u , 

(3) 



spectral functions corresponding to S , ^ TL (w, p) are found 



with projection operators P±(p) = (1 ± 707 ■ p)/2. The 
spectr; 
to be, 



--2ir\Z ± {p)8{u - uu±{p)) 
+ p\{u,p). 



h(p)s(u + ^(p))\ 



(4) 



In this limit, p I j ITL (u;,p) and p^ TL (uj,p) have two poles 
at uj = ±uj±(p) and ±u> T (p), respectively, with uj+(p) > 
and the residues Z±(p) — — p 2 )/(2m 2 ^) Q. Since 

the excitation spectra w±(p) becomes ttt-t at p = 0, 
tut is called the thermal mass. The spectral functions 
Eq. ([!]) have the continuum p±(u),p) in the space-like 
region, which originates from the Landau damping of 
the quark Q. Since the quark propagator Eq. |T]) is 
a function of u), p and tut, in the following we write 
5,^ L (oj,p) = S htl (lu,p) to show hit dependence ex- 
plicitly. 

While the HTL and one-loop approximations used in 
the above discussion may be invalid near T c , numerical re- 
sults from lattice QCD obtained large masses near T c Q. 
It is therefore plausible that the quarks have large ther- 
mal masses in the non-perturbative region near T c . The 
non-perturbative gauge interactions also produce an at- 
tractive interaction between quarks and antiquarks, and 
generate the chiral symmetry breaking. In order to model 
this possibility, we adopt Eq. fl} for the quark propaga- 
tor near T c and introduce a four-Fermi quark-antiquark 
interaction. We then arrive at the following model; 



(ipij 5 ip) 2 



(5) 



where ip is the quark field in the chiral limit, and Gs is 
the scalar coupling. The thermal mass m-r is introduced 
as a parameter to be varied by hand. We introduce the 
three-dimensional cutoff A to eliminate the ultraviolet 
divergence. This model with mj- = is equivalent to the 
Nambu-Jona-Lasinio model [ill. 



III. CHIRAL PHASE TRANSITION 

Let us first explore the chiral phase transition at finite 
temperature in the model of Eq. ([5]). In the mean-field 
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FIG. 1: Order parameters for the chiral phase transition a = 
-2G s (i>ip) for m T /A = 0, 0.1, 0.2 and 0.3 with fixed G S A 2 = 
13.01. The chiral transition is of second order irrespective of 
mr- For larger mr, o becomes smaller. 



approximation, the chiral condensate a = — 2Gs(ip' t P) is 
determined by minimizing the free energy 

rr 2 ^ f d 3 k 



4G S 



drk f duj 



(2tt) 3 



tanh - 



2T 



x trim In [S' R (w, fc; rriT, cr)] 1 , 
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where uj n = (2n + l)7rT is the Matsubara frequency for 
fermions and 



S(iu n ,k;m T ,a) = ([S™(iu n , k)}- 1 - a) 



(7) 



with S^, h being the Matsubara propagator of the 
quark in the HTL approximation; S^ L {iuj ni k) = 
S^ L (w, fe)| w _>j Wn . The retarded propagator S R in the 
right hand side of Eq. ^ is defined by the analytic con- 
tinuation S k (uj) = 5 , (it<j 9 T. ) ►c i j+ ir/ • The global mini- 
mum of Eq. ([(j]) satisfies the following stationary condi- 
tion (gap equation); 
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f d 3 k f dui oj p . 

J (27r)3 J ~ ^f trlmS (u,k;m T ,<T) 
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In general, Eq. ([5]) has several solutions corresponding to 
local extrema of Eq. ([6]). One must, thus, refer to Eq. (J6j) 
to find the global minimum. In our case, however, we 
have checked that Eq. ([8]) always has only one solution 
and it is certainly the minimum of Eq. ([6|). Therefore, 
we determine the order parameter only with Eq. (8). 

In Fig. [I] we show the chiral condensate a as a func- 
tion of T with several values of tot; although tut may 
change with a variation of T, we fix it to study the ef- 
fect of ttit- The scalar coupling GsA 2 = 13.01 is fixed 
so that T c /A = 0.3 at vht = 0. From the Figure, we see 
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that the chiral phase transition is of second order for any 
rriT- This result is natural in the sense of the symmetry 
argument of the phase transitions, since the symmetry 
of the system does not change with the incorporation of 
rriT- It is notable, however, that the order parameter a 
no longer has the meaning of the mass gap of the quark 
spectra with finite mf. 

One sees that, as mj- becomes larger, the chiral con- 
densate decreases and similarly the critical temperature 
does. This behavior is physically interpreted as the ef- 
fect of the excitation energy to^ in the quark spectra: Be- 
cause the chiral symmetry breaking is induced by the con- 
densate between the quark and anti-quark, the increase 
of the excitation energy of these particles suppresses the 
condensate. Here, we remark that if we incorporate the 
current quark mass m c in our model, the value of the 
chiral condensate becomes larger [ll|, although m c also 
causes a mass gap in the quark spectra. This is because 
m c explicitly breaks the chiral symmetry and works as 
an external field; in terms of the effective potential V(a), 
m c gives rise to a linear term proportional to m c a. 



IV. MESONIC EXCITATIONS 



mined by the analytic continuation II^-(u},p) 
n a {iv n ,p)\ iVn ^ u+iv with 

tr[S™(iu, m ,fc) 
x S^ h (iu n + iuj m ,p + k)] 



s,t=± 



d 3 k 

(27T) 3 



tr[P s (fc) 7 °P*(p + fc)7 ] 



dLL>idiL>2 f(ui) - f(u2) 
(27r) 2 iv n +UJ1-0J2 
x^ TL K,fc)p* TL (a, 2 ,p + fc) (10) 

and the Matsubara frequency for bosons v n = 2nirT. In 
Fig. [21 we show the diagrammatical representation of our 
approximation. 
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Since the chiral transition is of second order in our 
model, there should appear soft modes associated with 
the transition [3]. In this Section, we show that this is 
indeed the case. Using these soft modes, we then discuss 
the effect of the thermal mass on the mesonic excitations 
above T c . In the following, we limit our attention in the 
chirally symmetric phase above T c , where the excitation 
spectra in the scalar and pseudoscalar channels are de- 
generate. 

In the random phase approximation, the retarded me- 
son propagators in these channels D^(u>,p) is given by 



-1 



(9) 



where the polarization function ll^(ui,p) is deter- 



FIG. 2: Diagrammatic representation for the meson propa- 
gator D^(u>,p) and the polarization function li^(uj,p). The 
thick-solid line denotes the quark propagator in the HTL ap- 
proximation Eq. ([T]), while the thin-solid line represents the 
free quark propagator. 



The imaginary part of n^(u,p) is proportional to the 
difference between the decay and creation rates of the 
mesonic excitations. For p = 0, it is calculated to be 



ImZT*(w,0) = 
with 
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— (JppM + JpgM+IccM), (11) 
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— 7 + lf(u }+ )-f(cj + u; + )}p h _(u J + u }+ ,k) + 2ir J —-Z-[f(u>- - u) - - w, fe) 

— (bj <-> —u>), (13) 
1 f d 3 k r 

^ccM=- J dE[f{E) - f(u + E)]p\(E, k)p h _ {u + E,k)-(u>^ -w), (14) 



with f(E) — [exp(E /T) + 1] , f± = f(u>±), lo' ± — turn of normal quasiquarks or plasminos which satisfies 
du>±(p)/dp, and each p in Eq. (|12p denotes the momen- lu = — 2cj+, lu = 2lu- and lu = w~ — w+, respectively. 
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(a) 



FIG. 3: The diagrammatical interpretation for the decay 
processes included in Ipc(u) (left) and Zcc(^) (right). The 
dashed line denotes the mesonic excitations and the temporal 
direction is left to right. 
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FIG. 4: Each component of Imi7 R (w,p = 0) = Jpp(w) + 
/pc(w) + IccO*0 at T/A = m T /A = 0.3. 7pp(w), 7pc0*0 
and 7cc(^) include the decay processes into two quasiparticle 
poles, a quasiparticle pole and continuum, and two continu- 
ums, respectively. 



The first (second) term in Eq. (|12p includes the decay 
process into the normal quasi-quark and quasi-antiquark 
(the plasmino and anti-plasmino) , and the third term is 
the Landau damping of the mesons between the normal 
quasiparticle and the plasmino. The decay processes cor- 
responding to Ipc{to) and Ice 0*0 include the continuum 
/9 L (o;,p), i.e. the Landau damping of the quark. We show 
the diagrammatical interpretation of decay processes in- 
cluded in Ipc(w) and Ice 0*0 m Fig.OJa) and (b), respec- 
tively: Notice that the Landau damping is the scattering 
processes accompanied with a thermally excited particle. 

In Fig.[U we show each component of ImII R for a typ- 
ical parameter set T/A — ttit/A = 0.3. From the Figure, 
one sees that Ipp(w) has divergences at two energies ui\ 
and ui2, with wi/tot — 0.4 and u>2/mT — 1-8. These sin- 
gularities are due to the divergence of the density of states 
of the energy spectra 2w_ (p) and o>_ (p) — oj + (p) at finite 
momenta, and interpreted as the van Hove singularities 
EH 0] • While Ipp(cj) vanishes at u>\ < lu < ui2 due 
to the kinematics, IpcO*0 and Ice 0*0 take finite values 
in the whole range of energy. In particular, Ipc(aj) has a 
large contribution at to < 3mr. 

The excitation spectrum of the mesons is given by the 
spectral function 



■ (uj,p) = -2ImD R (cj,p) 



(15) 



In FigG3 we show p a (u>,p) for tot /A — Q and 0.3 at p = 
in the upper and lower panels, respectively; the temper- 
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FIG. 5: Spectral functions for hit = (upper) and mr/A = 
0.3 (lower) with several values of T above T c . The coupling 
Gs is determined so that T c /A = 0.3 for each mr- 



ature dependence of tot is not considered as before. In 
these Figures, the scalar coupling Gs is adjusted so that 
T c /A — 0.3 for each mx- In the upper panel, we see that 
there appears a peak in p a (uj,p), and the peak becomes 
sharper and moves toward the origin as T approaches 
T c from high temperature. This is the soft mode of the 
chiral phase transition Q ■ 

In the lower panel with finite %, p a (u>,p) behaves in 
a more complicated way. First, p^iuj^p) becomes zero at 
two energies lo\ and u>2 corresponding to van Hove sin- 
gularities in Ipp(w). Second, although there appears a 
sharp peak in p a (u>,p) as a soft mode of the chiral tran- 
sition, the heights of the peaks in p a (uj,p) are depressed 
compared from the upper panel with mx = 0. This re- 
sult contradicts the naive expectation that the mesonic 
excitation at u) < Itht becomes more tightly bound and 
thereby have a small width in this range. 

In order to understand the origin of the large decay 
width of the mesonic excitations with tot > 0, we show 
ImH CT (o;,p = 0) for tot/ A = and 0.3 in Fig. [6] One sees 
that the imaginary part with vht/A = 0.3 takes larger 
values at uj < 2tot than that with tot = 0, indeed. As 
discussed previously, the decay processes corresponding 
to IpcO*0 an d the van Hove singularities in Ipp(w) has 
a large contribution in this range. 



V. SUMMARY AND DISCUSSIONS 

In this work, we explored the effect of the thermal mass 
of the quark on the chiral transition and the mesonic 
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FIG. 6: Imaginary part of II n (uj,p = 0) at mr/A = and 
0.3 with T — T c and 1.4T C . After including mj, the imaginary 
part is enhanced at to < 2mr. 



excitations at finite temperature using a simple model. It 
is shown that the order of the chiral transition does not 
change by incorporating a thermal mass. We also found 
that the thermal mass reduces the chiral condensate. By 
studying the spectral function of the mesonic excitations 
in the scalar and pseudoscalar channels, we found that 
the peak of excitations in the spectral function is strongly 
suppressed below 2tot- We showed that the continuum 
in the quark spectrum, and the van Hove singularities, 
play a significant role to enhance the decay rate of the 
mesonic excitations below 2tot- 

The decay processes included in 7pc(w) and Iqq{io) 
contain the effect of the Landau damping of the quarks, 
i.e. the scattering process of the quasiparticles by a ther- 
mally excited particle as shown in Fig. [3] These processes 



are specific to nonzero T and it is natural that such pro- 
cesses are obtained in our calculation with 5^ L (w,p). 
Since the contribution of these terms, as well as the van 
Hove singularities, may not depend on the channels of the 
excitation qualitatively, our calculation strongly suggests 
that any mesonic excitations in the light quark sector are 
easily destroyed by the thermal effect in the QGP phase. 
It is notable that there can nevertheless appear sharp 
peaks in the mesonic spectra as the soft modes of the 
chiral transition near T c Q. 

Although we employed the fixed tot to show the nu- 
merical results, tot ~ gT may vary as a function of T. 
The energies u>\ and u>2, i.e. zeros of p a {u>,p), thus, 
change with respect to T. Moreover, the width of the 
collective excitation of the quark can depress the van 
Hove singularities. Therefore, the depression of the me- 
son spectra is not, unfortunately, an experimental signa- 
ture. 

While we employed the quark propagator Eq. ([1]) in 
this work, the quark spectrum near T c may have the large 
decay width and a complicated dispersion [l(| HE]. It 
is interesting to explore mesonic excitations with such 
quark spectra instead of . This is beyond the scope 
of this work. 
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